# Both lambda and delta move
max_W = 3000
W_values = seq(-max_W,max_W,by=10)
decalage = 0
lambdas = seq(4,2,length=20)
lambdas = 10**-lambdas
deltas = seq(4,2,length=20)
deltas = 10**-deltas
for (j in 1:length(lambdas)){
grass = sigmo_pluie(W_values,lambdas[j], decalage = decalage)
plot(
W_values,grass,
xlim=c(-max_W,max_W),
ylim = c(0,1),
main=paste("j=",j,"lambda=",round(lambdas[j],5))
)
abline(v=0,lty=2)
grass = sigmo_coeff_competition(W_values,deltas[10], decalage = decalage)
plot(
W_values,grass,
xlim=c(-max_W,max_W),
ylim = c(0,1),
main=paste("i=10",round(deltas[10],5))
)
abline(v=0,lty=2)
##############################################################################################
palette <- rainbow(length(deltas))
par(mfrow = c(1, 2))
plot.new()
legend("center", legend = paste("i=",1:length(deltas),"delta=",round(lambdas,4)),
col = palette, lty = 1,xpd = TRUE)
plot(
x=c(0),
y=c(0),
xlim=c(-max_W,max_W),
ylim = c(0,0.5),
xlab = "W", ylab = "",
main=paste("j=",j,"lambda=",round(lambdas[j],5))
)
for(i in 1:length(deltas)){
grass = sigmo_pluie(W_values,lambdas[j],decalage = decalage)*sigmo_coeff_competition(W_values,deltas[i],decalage = decalage)
lines(W_values, grass, col = palette[i])
abline(v=0,lty=2)
abline(h=1,lty=2)
abline(h=0.25,lty=2)
}
}








































